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We consider a coupled atmosphere-ocean model, which involves hydrodynamics, 
thermodynamics and nonautonomous interaction at the air-sea interface. First, we 
show that the coupled atmosphere-ocean system is stable under the external fluctua- 
tion in the atmospheric energy balance relation. Then, we estimate the atmospheric 
temperature feedback in terms of the freshwater flux, heat flux and the external fluctu- 
ation at the air-sea interface, as well as the earth's longwave radiation coefficient and 
the shortwave solar radiation profile. Finally, we prove that the coupled atmosphere- 
ocean system has time-periodic, quasiperiodic and almost periodic motions, whenever 
the external fluctuation in the atmospheric energy balance relation is time-periodic, 
quasiperiodic and almost periodic, respectively. 
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Abbreviated Title: Coupled Atmosphere-Ocean Dynamics 
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1 Introduction: A coupled atmosphere-ocean model 



The global ocean circulation consists of the wind-driven upper ocean circulation and a 
meridional overturning deep ocean circulation called the thermohaline circulation. The 
ocean thermohaline circulation involves water masses sinking at high latitudes and up- 
welling at lower latitudes. During the thermohaline circulation, water masses carry heat 
or cold around the globe. Thus, it is believed that the global ocean thermohaline circula- 
tion plays an important role in the climate |29j . 

The thermohaline circulation is maintained by water density contrasts in the ocean, which 
themselves are created by atmospheric forcing, namely, heat and freshwater exchange via 
evaporation and precipitation at the air-sea interface. Thus the thermohaline circulation is 
described by coupled atmosphere-ocean models (2711111 • Such coupled models also describe 
feedback of the thermohaline circulation on the atmospheric dynamics (e.g., temperature 
feedback) . 

Mathematical models are a key component of our understanding of climate and geophysical 
systems. The formulation and analysis of mathematical models is central to the progress of 
better understanding of the thermohaline circulation dynamics and its impact on climate 
change. 

We consider a coupled atmosphere-ocean model, with simplified atmospheric dynamics, 
i.e., the atmospheric dynamics is described by an energy balance model. 
This is a zonally averaged, coupled atmosphere-ocean model on the meridional, latitude- 
depth (y, z)-plane as used by various authors [SOI ESI IS |H1 HI El- This model has been 
shown to capture some interesting climate phenomena j3n | I35| H] . It is composed of a one- 
dimensional stochastic energy balance model proposed by North and Cahalan 25 , for the 
latitudinal atmosphere surface temperature 6{y, t), together with the Boussinesq equations 
for ocean dynamics in terms of stream function ip{y,z,t), and transport equations for 
the oceanic salinity S{y, z, t) and the oceanic temperature T{y, z, t) on the domain D = 



{(y,z):0<y,z<l}: 

et = eyy-{a + e) + Sa{y)-l{y)[So{y) + e-T]+ f{y,t),Q<y<l, (1.1) 



where q{y,z,t) = —Aip is the vorticity; velocity field is {v,w) = {ipz, —ipy)', a is a positive 
constant parameterizing the effect of the earth's longwave radiative cooling; Sa{y) and 
So{y) are empirical functions representing the latitudinal dependence of the shortwave 
solar radiation; j{y) is the latitudinal fraction of the earth covered by the ocean basin; Pr 



qt + J(g,V') 
Tt + J{T,tP) 
St + JiS,ij) 




AT, iy,z)GD 
AS, {y,z)eD. 





(1.3) 
(1.4) 
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is the Prandtl number and Ra is the Rayleigh number. The first equation is the energy 
balance model proposed by North and Cahalan |^. The forcing f{y, t) may arise from, for 
example, eddy transport fluctuation, stormy bursts of latent heat, and flickering cloudiness 
variables. Moreover, J{g, h) = Qxhy — gyh^ is the Jacobian operator and A = dyy + dzz is 
the Laplace operator. The effect of the rotation is parameterized in the magnitude of the 
viscosity and diffusivity terms as discussed in [55] . 

The no-flux boundary condition is taken for the atmosphere temperature 6{y,t) 

eyiO,t) = ey{l,t) = 0. (1.5) 

The fluid boundary condition is no normal flow and free-slip on the whole boundary 

-0 = 0, q = 0. (1.6) 

The flux boundary conditions are assumed for the ocean temperature T and salinity S. 
At top z = 1, the fluxes are specified as: 

Tz = So(,y) + {e-T)U=,, Sz = F{y), (1.7) 

with F{y) being the given freshwater flux. 
At bottom z = 0: 

= 5, = . (1.8) 

On the lateral boundary y = 0, 1: 

Ty = Sy = 0. (1.9) 

We also assume the following compatibility condition: 

S'M = S'oil) = F'{0) = F'{1) = 0. (1.10) 

The non-autonomous partial differential equation (jl.H) is only defined on the air-sea inter- 
face and it may be regarded as a dynamical boundary condition. The boundary condition 
()1.7() involves a coupling between the atmospheric and oceanic temperature at the air-sea 
interface. 

In the next section, we discuss the well-posedness of this coupled atmosphere-ocean model. 
Then we investigate the stability of this coupled system under external fluctuation in §3, 
atmospheric temperature feedback in §4, time-periodic, quasiperiodic and almost periodic 
coupled motion in §5, respectively. Finally, we summarize these results in the final section 
^6. 
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2 Mathematical Setup 



In order to use the standard result in (20| for the local existence, we homogenize inhomo- 
geneous boundary conditions for T, S on the top boundary z = 1 as in pH] . 
First, we construct two scalar functions: 

T: = f*r],{z), S: = S*r],iz), Ve G (0, 

where 

f* = [S,{y) + e]{l-e'-'), 

S* = Fiy)z, 

r]^{z) G C°°([0, 1]) is given by 



1, l-e< z <1, 

rie{z) = ^ increasing, 1 — 2e<z<l — e, 
0, < z < 1 - 2e. 



Then, set 



T = T-T* S = S-S* 



By ()1.5() and ()1.1U|) . we see that the boundary conditions (jl.Tf) for the new variables T and 
S are homogenized and do not affect other boundary conditions. Thus (|l.lj) - (|1.9jl become 
(for the simplicity, we still use T and S instead of T and S) 

et = eyy-{a + e) + Sa{y)-j{y)[So{y)+9-T] + f{y,t), < y < 1, (2.11) 
qt + J{q,i')=PrAq + PrRa{Ty-Sy + T:y-S:y), {y,z)eD, (2.12) 

rt + j(r,^) + j(t;,v) = Ar+ 

[((1 - e'-^il - 7(y)) - e'^n^.iz) + (1 - e'~^W:{z) + 2e'~^r^i{z) - e'-^r^,]e 

-(l-ei-^)r/,(z)7(y)r(y,l)+5, {y,z)GD, (2.13) 
St + J{S,i;) + J{S:,tP)=AS 
+ F"{y)zrjAz) + F{y){2rj'Jz) + Z7]'^{z)), {y, z) G D, (2.14) 

where 

g{y, z,t) = -{l- ei--)r?,(z){-a + 5,(y) - -i{y)SM + f{y, t) - S'^iy)} 
+ [(1 - e'-nri':iz) + 2e'-Xiz) - e'-^r^,]So{y). 
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The corresponding boundary conditions become: 
On the whole boundary, the fluid flow satisfies 



^ = 0, q = 0. (2.15) 

The boundary conditions for the atmosphere temperature 9{y, t) (defined only on the 
air-sea interface) are 

eyiO,t)=ey{l,t)=0. (2.16) 

The boundary conditions for the ocean temperature T and salinity S become 
At top z = 1: 

r, + r|,=i = 0; S, = 0. (2.17) 

At bottom z = 0: 

T, = S, = 0. (2.18) 

At the lateral boundary y = 0, 1: 

Ty = Sy = 0. (2.19) 
The appropriate initial data Oq^qq^Tq, Sq are also assumed. 

Using the theory in [20], we can obtain the following local existence theorem for problem 
(2.1)-(2.9) (that is (1.1)-(1.9)). 

Theorem 2.1 (Local Well-Posednessj Let Oq € H'^{0, 1), go S H^{D), Tq, Sq G H^iD), 
f G L°°(0, (X); L^(0, 1)). Assume that the physical data satisfy 7(7/) G L°°(0, 1), and 
So{y),F{y) G H^{0,1) and aslo assume that the compatibility condition U.l(]\) be sat- 
isfied. Then the coupled atmosphere- ocean system \2.11\) ~ Ji2.iy\) (that is hl.l\) - \l.y\) ) has 
a unique (The uniqueness of S is up to a constant) local solution satisfying 

e G L~(0, r; 7/1(0, 1)) nL2(0,r;F2(0,l)), 

q G L^{0,T;H^{D))nL\0,T;H\D)nH}^{D)), 
T, 5 G L°°(0, r; H^{D) x H^{D)) f] L^{0, r; H^{D) x H^{D)), 
where r depends on initial data and physical data Sa{y), So{y), F{y) and f{y,t). 

Since —Aip = g G Hq, we get ip G Hq{D)i^H^{D). Hence the Jacobian J(-, •) is continuous 
from H^{D) x H^{D) L^{D) x L^{D). 

In order to obtain the global existence, we need a priori estimates. First, we give a priori 
estimates for (|1.1|) - (|1.9|) in L^. In the sequel, || • || and || • ||i denote the norm of and 
respectively. 
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Multiplying by 6 and performing the integration by parts, we conclude that 

~\\ef = -\\0yf-\\of 

- a 9dy + r SaOdy - r{y) [So + O- T]9dy + [' fOdy, (2.20) 
Jo Jo Jo Jo 

By the Cauchy-Schwarz inequality, we arrive at 



1^ 
2Jt 



< -\\Oyf + Ci\\0f + 6i|lT(y, l)f + Ml, (2.21) 



and 



Note that 



where constant Mi depends on HS'all, ||S'o||, a and sup ||/||, constant Ci depends on 

0<t<oo 

||7||loo and ei > 0, ei > will be chosen later. 

Multiplying (|1.2() by q, performing the integration by parts and using the property of 
Jacobian and H1.6|l . we have 

~hf = -PrW^qf + PrRa jjTy - Sy)q. (2.22) 
Similarly, from (|2.13)) and (|2.14|) . we have 

~\\Tf = -llVrf + j\So + e - T{y, l)]T{y, l)dy, (2.23) 

-\\VSf+ [' Fiy)S{y,l)dy. (2.24) 
Jo 

PrRa j^iTy - Sy)q < + ^^^^{\\Tyf + \\Syf) 

< ^\\Vqf + PrRa^XiiWVTf + \\VSf), 

where Ai is a constant in the inequality < Ai||Vw||,w G Hq. Thus (|2.22|) can be 

rewritten as 

2jM - "T"^'^"' + PrRa'XiiWVTf + \\VSf). (2.25) 
Multiplying (fT^ by 2PrRa?\i and fTM by 2PrRa^\i, and adding to flT^ . we get 

]^^^{\\qf +2PrRa'\i{\\Tf + \\Sf)) 
Pr 

< -^ilV^f - PrRa^XiiWVTf + \\VSf) 

+ 2PrRa^Xi f' [So + O- T{y, l)]T(y, l)dy + 2PrRa^\i f' F{y)S{y, l)dy. (2.26) 
Jo Jo 
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By the Cauchy-Schwarz inequality and the trace ineauaUtvf|13jl. we have 

-^{\\qf + 2PrRa^Xi\\Tf +2PrRaHi\\Sf) 
2 at 

Pr 

< -^llVgf - PrRa^\i{\\VTf + \\V Sf) 

+ C2\\ef - ^||T(y, l)f + e2{\\VSf + \\Sf) + M2, (2.27) 

where C2 depends on Pr, Ra and Ai, and M2 depends on Pr, Ra, Ai, ||5o|| and ||F||. Choos- 
ing ei < i and €2 < ElM^, combining f^Tm with ^^TITl . we obtain 

^■^iPf + \\qf + 2PrRa'^Xi\\Tf + 2PrRa^\i\\Sf) 

< -a{\\vef + \\vqf + iivrf + \\vsf) + C3{\\ef + \\sf) + M3, (2.28) 

where C3 depends on Ci and C2, M3 depends on Mi and M2, and a is a positive constant 
depending on Pr, Ra and Ai. By the Gronwall inequahty, we have 

\\0f + \\qf + \\Tf + \\Sf 

rb 

+ (llV^f + ||Vgf + ||VT||2 + ||V5f)dt<C7i(6), (2.29) 
Jo 

for any given future time 6(0 < b < 00, < t < b) and some positive constant Ci{b) 
depending on 6, C3 and M3. By a similar argument in J3] (here we need to obtain the 
estimates in x x x x for system of (|2.11|) - (|2.14|) in order to avoid the 
trouble of non-homogeneous boundary conditions, we omit the details here since we will 
give a similar proof in §4), we have 

\\vef + \\vqf + iivrf + ||V5f + 

[\\\Aef + \\Aqf + \\ATf + \\ASf)dt < C2{b), (2.30) 
Jo 

for any given 6(0 < 6 < 00) and some positive constant C2(6) depending on j|5o||/f2, 
\\F\\h2, 6 and Ci(6). 

Remark 2.2 In fact, the estimate we get in \2. is for T and S, from which we can 
get the estimate for oringinal T and S using the estimate for 9. 

Remark 2.3 Since the equivalence ||r|p-|-|| ATp with ||T|||^2 and the same for S , together 
with irm . we can replace f^{\\Aef + \\Aqf + \\ATf + \\ASf)dt hy f^{\\e\\]j2 + \\q\\]j2 + 
\\T\?H2 + \\S\W2)dt in the estimate JO?^) . 
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With these global estimates, we have the following global existence theorem for the coupled 
atmosphere-ocean system: 

Theorem 2.4 ("Global Well-Posednessj Let Oq G H^{Q,l), Qo e H^{D), To,So G 
H^iD), f G L°°(0,oo;L2(0, 1)). Assume that the physical data satisfy 7(y) G L°°(0, 1), 
and So{y), F{y) G H^{0,1) and aslo assume that the compatibility condition \1.1U\) be 
satisfied. Then for any given 6(0 < b < oo), the coupled atmosphere- ocean system h2.11]) - 
\2.iy\) (that is hl.l]) - \l.y\) ) has a unique (The uniqueness of S is up to a constant) global 
solution satisfying 

e G L°°(0, 6; i/i(0, 1)) n l2(0, 6; //^(0, 1)), 
q G L°°(0, 6; HI{D)) n L\Q, b; H\D) D i^o'P)), 
T,S & L^iO,b;H\D) x H\D)) D L'^iO,b; H^{D) x H^{D)). 

In the rest of this paper, we assume the conditions in this theorem are satisfied, so that 
we always have global unique solutions. 

In the next section, we consider the stability of the above coupled atmosphere-ocean 
system with respect to the external fluctuation f{y,t) in the atmospheric energy balance 
dynamics 

3 Stability under External Fluctuation 

Paleo-evidence on the instability of the thermohaline circulation is now abundant. Numer- 
ical work suggested that a sufficiently large external forcing (such as external fluctuations 
in the atmospheric energy balance model and the freshwater flux at the air-sea interface) 
could destabilize or shutdown the thermohaline circulation j2H]. This indicates that cur- 
rent capacity of carrying heat poleward by the thermohaline circulation may change when 
the freshwater budget is altered. Since the thermohaline circulation's important role in 
redistributing the heat around the globe, a breakdown or instability of the current ther- 
mohaline circulation may lead to dramatic climate change j^H]- Because of this relation 
between the thermohaline circulation and climate change, there is growing interest in its 
stability or instability. This makes the stability issue of the thermohaline circulation not 
only of scientific but also of great practical importance. 

In this section, we prove the stability of the coupled atmosphere-ocean system with respect 
to the external fluctuation f{y,t) in the atmospheric energy balance dynamics i.e., 
the continuous dependence of solution on f{y,t) in the space H^. 
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Assume that {6i,qi,Ti, Si} and {92,q2,T2, S2} are solutions with respect to fi{y,t) and 
f2{y,t). Let 

9 = 01-02, q = qi-q2, f = Ti- T2, S = Si- S2, / = /i - /2, 

then 0, q, T and S satisfy 

Ot = Oyy-0 - j{y)[0-f]+f{y,t), 0<z<l, (3.1) 

qt + J{qi, V'l) - J{q2, V'2) = PrAq + Pri?a(2; - Sy) , {y, z) G L>, (3.2) 

2i + J(ri,Vi)- J(r2,V2) = AT, (y,z)GA (3.3) 

5t + J(5i,V'i)- J(52,V'2) = A5, (y,z)GD, (3.4) 

The corresponding boundary conditions are as follows. 
On the whole boundary: 

^ = 0, AV' = g = 0. (3.5) 



0y{^,t) = ~0y{\,t)=^. (3.6) 

At top z = \: 

f, + f\,=i = 0; S, = 0. (3.7) 

At bottom z = 0: 

f, = S, = 0. (3.8) 

At the lateral boundary y = 0,1: 

fy=Sy = 0. (3.9) 

Similar to the discussion in §2 above, we have 

^|ll^T = -ll^".f-II^T 

- /' r{y)\0\^dy + [' ^{y)T{y, l)0dy + [' f0dy, (3.10) 
Jo Jo Jo 

IjjQf + jj.J{qi,^i) - Jiq2,^2m = -Pr\\Vq\\^ + PrRajjJy - Sy)q, (3.11) 
\jt\\n' + JjJ{Ti,A)-J{T2,^2))f = -\\Vff + j\0-ny,l)]f{y,l)dy, (3.12) 

~\\~S\? + J^{J{Si,^i) - J{S2,^2))S = -\\VSf. (3.13) 

In order to estimate the terms Jj^{J{qi,ipi) — J{q2,'tp2))q, loi'^^'^^^'^'^) ~ •^(^2, '^2))T and 
Jj-,(J(Si,jpi) — J{S2,ip2))S, we need the following lemma: 
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Lemma 3.1 The nonlinear Jacobian operaror J{u,v) has the following property 

\\J{UI,U2) - J{VI,V2))\\ < 

(llVuill + \\Vn2\\ + ||Vt;i|| + ||Vr;2||)(||V(ni - ^;i)|| + ||V(n2 - v^W, (3.14) 
for every Ui,Vi € H^{D) (i = 1, 2). 

The proof of this lemma is in . 
By Lemma 3.1, we have 

/ (J(gi,^i) - J((?2,V'2))g < (1 + Ai)2(||Vgi|| + 

JD 

f (J(ri,v^i)- j(r2,^2))T 

JD 

< (llVTill + IIVT2II + AidlVgill + ||Vg2||))(||Vr 

and 

/ (J(5i,Vl)- J(52,V'2))5 
JD 

< (||V5i|| + IIV52II + AidlVgill + ||Vg2||))(||V5 

Note that 

PrRa [ {fy - Sy)q < ^\\qf + PrRa^\i{\\fyf + \\Syf) 

< ^llVgf + PrRa^XiiWVff + \\VSf). 
Then by a similar argument as in §2 and using the Cauchy-Schwarz inequality, we conclude 
that 

+ Uf + 2PrRa^\i\\ff + 2PrRa^Xi\\Sf) 

^ (XL 

< CiiWef + \\qf + 2PrRa^Xi\\ff + 2PrRa^Xi\\Sf) + \\ff, (3.15) 

where C4 depends on Pr,Ra, Ai, ||7|[l°° as well as the H^—novm of q,T and S. By the 
Gronwall inequality, we further have 

\\9f + \\qf + ||rf + \\Sf < C{b)\\f\\\ (3.16) 

for any given b (0 < b < 00, < t <b) and some positive constant C{b) depending on b and 
C4. Furthermore, we can obtain the similar estimates for the gradient of {9,q,T, S}, we 
omit the proof here, as the similar derivation will be done in §5. Thus the solution differ- 
ences 6, q, T, S and / are bounded when the external fluctuation difference / is bounded. 
So we have the following stability theorem. 

Theorem 3.2 (^Stability under the external fluctuation^ The coupled atmosphere- 
ocean system h2.11\) - \2.1Ji\) (that is il.l\) - (Tm )) is stable under the external fluctuation in 
the atmospheric energy balance model. Namely, the solution of the coupled system depend 
continuously on the external fluctuation f in . 



Vg2||)||Vg||||g| 

+ Ai||Vg-||)||r 

+ Ai||Vg-||)||5 
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4 Dissipativity and Atmospheric Temperature Feedback 



The ocean and the atmosphere are constantly interacting through the air-sea exchange 
process. It is expected that the thermohahne circulation could provide feedback to the air 
temperature. This is a direct impact of the thermohahne circulation on the climate. It is 
desirable to predict or estimate this feedback. 

To this end, let us estimate the air temperature 9 in the mean-square norm, in terms 
of the freshwater flux F{y), external fluctuation f{y,t) in the energy balance model, the 
earth's longwave radiative cooling coefficient a, and the empirical functions Sa{y) and 
So{y) representing the latitudinal dependence of the shortwave solar radiation, as well as 
physical parameters Pr and Ra. 

We will also show that the system generated by ()1.1|) - H1.9|) is a dissipative system in the 
sense of ^2] or under some conditions, that is all solutions {9, q, T, S} enter a bounded 
set (so-called absorbing set) in H^{0,1) x Hq{D) x H^{D) x H^{D) after a finite time. 
Since 

— / Sdydz = / F{y)dy = constant. 
dt Jq Jo 

For simplicity, we assume that 

= 0, / Sdydz = {) (4.1) 
JD 

and 

< 7(y) < 1 or < 7(y) < 1. (4.2) 

First, we derive a uniform estimate for {9,q,T,S} in L2(0, 1) x L'^{D) x ^^(D) x L'^{D). 
Using the standard energy estimate as given in §2, we get 



l_d 

2dt' 



9 IP 



, f\dy+ 1^ Sa9dy- f\iy)[So + 9-T{y,l)]9dy+ f' f9dy 
Jo Jo Jo Jo 

< -\\9yf - {I - e)\\9f - M My)\\9f 

y6[o,i] 



+ -[a^ + \\Sof + \\Saf + Wff] + \\riy)\\Loo [' \9\\T{y, l)\dy. (4.3) 
e JO 



^^^{\\qf +2PrRa'X,{\\Tf +\\Sf)) < 
-^llVgf - PrPa^AidlVrf + \\VSf) 
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+2PrRa''Xi [\So + 9 - T{y,l)]Tiy,l)dy + 2PrRa''Xi (\{y)S{y,l)dy. 
Jo Jo 

< -^llVgf - PrPa^XiiWVTf + \\\/Sf) - 2PrRa^Xi{l - e) /' \T{y, l)\^dy 
2 Jo 

+ 2PrRa'X, f \e\\T{y, l)\dy + ^i:^^\\S^f + P^'^^'^^l \\Ff ^ ^^^^s{y, l)f, (4.4) 
Jo e ei 

here e > and ei > will be chosen later. By the trace inequality, we have 

\\Siy, l)f < C(||V5f + \\Sf) < C{1 + Ai)||V5f , 
where Ai is the constant in the following Poincare inequality (note that J^, Sdydz = 0) 

<Ai||V5f. (4.5) 
Choosing ei = 2(1+11)0 ' ^^^^ 1|4.4() can be written as 

^±{\\qf +2PrRa'Xi{\\Tf +\\Sf )) < 

-^llVdP - PrPa'XiiWVTf + ^||V5f ) - 2PrRa''Xi{l - e) J^' |r(y, l)fdy 

+ 2PrRa^Xi \e\\T{y, l)\dy + ^""^"'^^ + 2PrRa^Xiil + Xi)C\\Ff. (4.6) 

JO e 

Then, multiplying (|4.3() by 2PrRa^Xi and adding to H4.6p . we get 

|2 + 2Pri2a2Ai(||^f + llTf + ||5f )) < 

Pr 1 

-2PrRa'Xi\\eyf - —\\Vqf - PrPa^Aidl VTf + i^W'^Sf) 

-2PrRa^X^{{l-e)\\ef- ini My)\\Of-{l+h\\L^) [' \T{y,l)\\e\dy+{l-e) [' \T{y,lfdy) 

ye[o,i] Jo Jo 

+ :^!j^[a2 + + +\\Saf + Wff] + 2PrRa'X,{l + X^)C\\Ff. (4.7) 

By (|4.2|) . when < 7 < 1, we could choose e such that (since inf 7(7/) = now) 

j/e[o,i] 

4(1 - ef > (1 + ||7|U^)^ i.e. 6 < := oq. 



For example, we choose e = then 

-2PrRa^Xii{l - e)\\0\\ - (1 + hb-) /' \T{y, l)\\9\dy + (1 - e) /' \T{y, l)\^dy) 

Jo Jo 
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<-^^^^^m? + \\nyA)f). 

If < 7(y) < 1 as in (|4.2|1 . we denote inf ^ 7(7/) = (3q. Then we take e = ^ to obtain 

-2PrRa''Xi{{l - e + /3o)||0|| - (1 + UtIIl-) /' |T(y, + (1 - e) /' |r(y, l)\''dy) 

Jo Jo 

D 

So, in the case of < 7(7/) < 1, (|4.7p can be written as 

^ 2 + 2Pri?a2Ai(||ef + ||rf + \\Sf)) < 



2dt 

Pr 1 

-2PrRa^X,\\eyf - —\\Vqf - PrRa^Xi{\\VTf + -||V5f ) 

+ ^!:^[a2 + 5 |5^||2 + |[5^||2 + ||^||2j ^ 2Pri?a2Ai(l + X^)C\\F\\\ (4.8) 
ao 4 

For < 7(y) < 1, we wiU have similar estimate. Since 

T\y,z)-T\y,l) = 2jyT,dz, 

we further have 

llTf <2 f\Tiy,l)\^dy + i\\\7Tf. (4.9) 
Jo 

Using the Poincare inequhty again for q and letting ai = min{-^^, |j ^} and /3 = 
uim{^, ErRn^h»a^ PrRfx, }^ the estimate (4.8) becomes 

^^^{\\qf + 2PrRa'X,{\\ef + \\Tf + \\Sf)) < 

-a^i\\qf + 2PrRa^X^{\\9f + \\Tf + \\Sf) 
-mOyf + llVgf + llVTf + \\VSf + |jT(y, l)f) 

+ :^!j^[a2 + ^||5of + ll/f ] + 2PrRa'Xi{l + Ai)C||Ff. (4.10) 

Using the Gronwall inequality, we finally obtain the mean-square norm estimate for the 
solution of the coupled atmosphere-ocean model: 

\\qf + 2PrRa^X,m\'+\\Tf + \\Sf)< 
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e-"^*(||gof + 2PrRa^\i{\\eo\? + W^f + \\Sof)) 

+ ^^^^[-(a^ + hSoW + \\Sa\? + ll/f ) + 2(1 + Ai)C||Ff ] (4.11) 
ai ao 4 

In particular, we get the mean-square norm estimate for the atmospheric temperature 
feedback 

zai ao 4 

This atmospheric temperature feedback estimate is in terms of physical quantities such 
as the freshwater flux F{y), external fluctuation f{y,t) in the energy balance model, the 
earth's longwave radiative cooling coefficient a, and the empirical functions Sa{y) and 
So{y) representing the latitudinal dependence of the shortwave solar radiation, as well 
as physical parameters Pr and Ra. Here Ai and Ai are the constants in the Poincare 
inequality on the domain D in the cases of zero Dirichlet boundary condition and zero 
mean value, respectively. Moreover, C is a constant depending only on the domain D, 
"0 = «i = min{i^, f, i, ^}. 

We can furthermore derive solution estimate in norm. To do so, we first get from 

f'^^W^yW^ + ll^^f + llVrf + IIV^IP + \\T{y, l)|p) < 
-^e-»*(||(?of + 2PrRa^X,{\\eof + \\Tof + \\Sof)) 

+ ^^'^f^^ (-(«^ + f + \\Saf + Wff) + 2(1 + XM\Ff). (4.13) 

So, let ||(7o|p + 2Pri?a^Ai(||6'o|p + llTolP + Ij^olP) be bounded by some (big) upper bound 
i?2 and denote W = ^^^^(^-(a^ + |||5of + \\Saf + Wff) + 2(1 + Xi)C\\Ff). Then 
there is a time t* > — In such that 

— ai M 

\\qf + 2PrRa^Xii\\ef + ||rf + \\Sf) < 2M^, t > t* (4.14) 

and 

r\pyf + llVgf + llVrf + \\VSf + \\T{y, l)f ) < J-M^, t > t*. (4.15) 
Jt ■ Pi 

Next, we derive a uniform estimate of gradient of {6, q, T, S} in L^(0, 1) x (D) x L^{D) x 
LP'{D). In order to avoid the difficulty caused by the non- homogeneous boundary condi- 
tions, we use equations (|m]) ~ (|TT^ instead of (P~T|) - (fO)) . 



15 



Multiplying (|2.11|) ~ H2.14|) by —Oyy, — Ag, —AT and —AS" respectively, integrating over 
(0, 1) and D, noting that S* G H^{D) is known and T* is only dependent on 6, we get 

+ a f' dyydy - [' SaOyydy + f' r{y)[So + 9- T{y, l)]9yydy - f' fOyydy, (4.16) 
JO Jo Jo Jo 

~\\Vqf = -Pr\\Aqf + J^J{q,^P)Aq- J^PrRa{Ty - Sy + - S;)Aq, (4.17) 

i^dlVrf + \\Tiy, 1)||2) = -\\ATf + 1^ J{T + T\^)AT 

- [ [((1 - e'-^l - 7(y)) - e'-^)r^,{z) + (1 - e'-^W^iz) + 2e'-^r^Uz) - e'-^r^,]eAT 
Jd 

+ f {l-e^-')7],{z)j{y)T{y,l)AT- f gAT, (4.18) 
Jd Jd 

^|l|V5f = -||A5f + J{S + S*, ^)AS 

- [ F"{y)zr^,{z)AS- ( F{y){2r^[ + ztj," {z))AS. (4.19) 
Jd Jd 



Note that 



a I Oyydy - I SaOyydy+ [ r{y)[So + 6 -T{y,l)]9yydy - [ f9yydy 
Jo Jo Jo Jo 

< jPyyf + + \\Saf + \\Sof + Wff + \\9f + Wny, l)f ], (4.20) 
- PrRa{Ty - Sy + T; - S*y)Aq < 

^\\Aqf + ^^(llVTf + ||V5f + \\9yf) + ^-^^{\Kf + WF'f), (4.21) 
[((1 - e'-^l - 7(y)) - e'-^)r^,{z) + (1 - e'-^Veiz) + 2e'-^7^[{z) - e'-%]9AT 

D 

+ f {l-e'-')ve{zh{y)T{y,l)AT- f gAT 
Jd Jd 

< l||Ar||2 + C(a2 + \\Saf + \\So\\j,2 + ll/f + \\9f + \\Tiy, l)f). (4.22) 

- / F"{y)zr^,{z)AS- f F{y){2r^', + zr^J' {z))AS < ^ASf + C\\F\\l,. (4.23) 
Jd Jd 2 

About the estimates of J{q, Tp)Aq, Jj^^ J{T + T*,^)AT and J{S + S* , ^j)AS, similar 

to JH], we have 

' J(g,^)A<7 < CIIAV'IIIIVqIIIIAc^II = C||g|| ||Vg|| ||Ag||, (4.24) 

D 
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/ J{T + T*,iP)AT < C||g||(||AT|| + ||Vr|| + ||r||)||Vr|| + ||g||(||0j^|| + ||5;||)||Ar||, (4.25) 
Jd 

[ J{S + S*,il^)AS < C||g||(||A5|| + \\VS\\ + ||5||)||V5|| + ||g||||F'||||A5||. (4.26) 
Jd 

Using the Cauchy-Schwarz inequality, H4.2Up - (|4.26|) and (|4.14|1 . when t >t*, we get 
IftiW^W + llV'zf + SPrRa^XiiWVTf + \\T{y, l)f + ||V5f )) 

< Cipyf + \\Vqf + 3PrRa^Xi{\\VTf + \\T{y, l)f + \\VSf) 

+ C{a^ + \\Saf + \\So\\l2 + Wff + \\F\\l,). (4.27) 
By (|4.15|) and a uniform Gronwah lemma (j^-ilj). we obtain 

\\Gyf + llVgf + 3PrRa^Xi{\\\7Tf + \\T{y, l)f + \\VSf) 

< C{a^ + \\Saf + \\SofH2 + Wff + \\F\\l-2). (4.28) 

By (14.141) and 1)4. 28() . we know there exists an absorbing sets B in H^{0, 1) x Hq{D) x 
H^{D) X H^{D) for the solution of irm)- (fT"9|l : 

B={{9,q,T,S} : \\e\\l+\\q\\l+\\T\\l + \\S\\l < C{a' +\\Saf + \\SofH2 + \\ff +\\F\\l,)}, 

(4.29) 

that is, for every bounded set in H'^{0, 1) x H^iD) x H^{D) x H^iD), when t > t* + 1, 

the solution of Hl.l() - ()1.9() will enter into the B. 

We summarize our results in section in the following theorem. 

Theorem 4.1 ^Atmospheric temperature feedback and dissipativityj Assume 
that the freshwater flux F{y) has zero mean as in J^. j| ) and the ocean basin's latitudinal 
fraction function ■j{y) is bounded as in i4.S\ ). Then the coupled atmosphere- ocean system 
m]\)-(r^ has an absorbing set in (0, 1) x (D) x H'^ (D) x H'^ (D) as given in 
More importantly, the atmospheric temperature feedback 0(y,t) is bounded in mean-square 
norm in terms of physical quantities such as the freshwater flux F{y), external fluctuation 
f{y,t) in the energy balance model, the earth's longwave radiative cooling coefficient a, 
and the empirical functions Sa{y) and So{y) representing the latitudinal dependence of the 
shortwave solar radiation, as well as the Prandtl number Pr and the Rayleigh number Ra 
as in i4.1S\ ). 

Remark 4.2 Due to J^-, S = 0, and J^, So = 0, we obtain J^S = 0. So, as seen in 
the discussion of this section, we know that the stability (proved in §5j under the exernal 
fluctuation is uniform in time t when < "y{y) < 1 or < 7(2/) < 1. 



17 



5 Strong Contraction and Almost Periodic Atmosphere-Ocean 
Dynamics 



In this section, we study the coupled atmosphere-ocean dynamical response to almost pe- 
riodic (in particular, periodic and quasi-periodic) external fluctuation f{y,t) in the atmo- 
spheric energy balance model A central question is: Does the coupled atmosphere- 
ocean system respond almost periodically to almost periodic external fluctuation f{y,t)? 
To answer this question, we need to understand the strong contraction property of the cou- 
pled atmosphere-ocean system in the absorbing set B defined in (|4.29|) . Let {9^,q'^,T'^, S*} 
be two trajectories corresponding to initial values {6Q,qQ,T(j, Sq} G S for z = 1,2. Note 
that these trajectories remain inside as i? is a forward invariant set. Their difference 

56 = 9^ - 9^, 5q = q^ - 5T = T^ - T^, 5S = - 

satisfy the following equations: 

69t = 69yy-69 - -f{y)[69 - 6T], < z < 1, (5.1) 

dqt + J{qi,^i)-J{q2,^2) = PrA6q + PrRa{6Ty-6Sy),{y,z)eD, (5.2) 

5rt + J(ri,Vi)- J(T2,V2) = A6T,{y,z)eD, (5.3) 

6St + J{Su^l)-J{S2,^P2) = ASS,iy,z)eD, (5.4) 

The corresponding boundary conditions are: 

5ey{0,t) = 69y{l,t) =0. (5.5) 



On the whole boundary: 
At top z = 1: 
At bottom z = 0: 
At the lateral boundary y = 0, 1: 



5^jJ = 0^ Sq = 0. (5.6) 
6n + STU=i = 69; 6S, = 0. (5.7) 
6T, = 6S, = 0. (5.8) 



6Ty = 6Sy = 0. (5.9) 

The initial conditions are: 

69o = 9l- 9l 5qo = ql- ql 5To = 4 - Tl 6So = - Sl (5.10) 
where 59q = 59{y, 0), 6qQ = 5q{y, z, 0), 5Tq = 5T{y, z, 0) and 5Sq = 6S{y, z, 0). 
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Using energy estimates as in §3, we have 

^ \\5ef = -\\5eyf - \\69f - l^j{y)m'dy - j^^{y)5e5T{y, l)dy, (5.11) 



1 d 
2di 



2dt 

2 



-Pr\\Vdqf-J {J{q^,%b^)-J{q'^,ip^))6q- J PrRa{6Ty - 6Sy)dq, (5.12) 



~\\STf = -\\V5Tf 



+ [\5e-5T{y,l)my,l)dy- [ {J{T\'^^) - J{T\'^^))5T, (5.13) 

JO JD 

^ = -||V55f- / {J[S\i^^)-J{S\i^^))5S. (5.14) 

JD 



D 



2dt 

Using Lemma 13.11 we imply 

{Jiq\i^') - Jiq^i^'))6q < \\Jiq\i^') - Aq' ,^l^^)\\\M < 

D 

\Vq'\\ + \\S/q^\\ + \\V^P'\\ + || V^2||)^||^^^|| ^ || V.J^'ll) ll-^gll 
<{l + X,f^,{\\q'\\ + \\q'\\)\\V6qf, 
{J{T\^') - J{T^^'))6T < \\J{T\i;') - J{T\^')\\\\6T\\ 

< (||vri|| + iivT^ii + ||v^i|| + ||vv^2||)(|[^^^|| ^ w^smm 

< (||VTl|| + IIVT^II + X,i\\Vq'\\ + WVq'WMWVdTW + X,\\V6q\\)\\5Tl 

- / {J{S\i;') - J{S^i;^))6S < \\J{S\i;') - J{S\^')\\\\6S\\ 

JD 

< iwvs^w + iivs^ii + iiv^^ii + \\vi>^\\){\\v5S\\ + IIVJV'IDII^^II 

< {\\VS'\\ + WVS^W + \i{\\Vq'\\ + ||Vg2||))((|jv55|| + Ai||V5g||)||55||. 

Other terms in (|5.1iP " H5.14p can be estimated as in the proof of the existence of the 
absorbing set in the last section. So we have 

~m\? + 2PTRa'\,{\\5ef + \\5Tf + \\5Sf)) < 



aiiWSqf + 2PrRa'\i{\\5er + H^Tf + WSSW')) 
'^0yf + l|V5gf + \\V6Tf + \\\V5Sf + ||5r(y, ] 
+ {l + X,f^^{\\q^\\ + \\q^)\\V5qt + 



2PrRa^\i{\\VT^\\ + \\VT^\\ + \i{\\Vq^\\ + || Vg2||))((||VJr|| + Ai || V5(/||) ||5r|| 
+2PrRa^\i{\\VS^\\ + \\VS^\\ + \i{\\Vq^\\ + \\V q^\\)){{\\V 5S\\ + Ai||VJ(7||)||<^5i|. 
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Now we assume that C(a^ + \\Sa\\'^ + 115*011^2 + + ||-^|||^2) is small enough. This is 
a condition imposed on the physical quantities such as the freshwater flux F{y), external 
fluctuation f{y,t) in the energy balance model, the earth's longwave radiative cooling 
coeflicient a, and the empirical functions Sa{y) and So{y) representing the latitudinal 
dependence of the shortwave solar radiation. 
By the Cauchy-Schwarz inequality, ()4.5() . 1)4. 9() and 1)4. 29() . we get 

~{¥qf + 2PrRa^X,{\\5ef + \\6Tf + \\6Sf)) < 

-ai{\\6qf + 2PrRa'^Xi{\\5ef + \\6Tf + \\6S\\^)) 
- Y^W^Oyf + \\V5qf + \\V6Tf + ^\\W5Sf + \my, l)f). (5.15) 

By the Gronwall's inequality, we obtain the strong contraction in L^(0, 1) x L^{D) x 
LP'{D) X L^{D) for the solution of the coupled atmosphere-ocean system (|1.1|1 - (|1.9|1 . that 
is 

\\Sqf + 2PrRa^M{\\def + \\6Tf + \\5Sf) < 
e-°^*(||<5gof + 2PrRa''X,{\\5eof + \mf + \\6Sof)). (5.16) 

Next, we can show the strong contraction of gradient in L^(0, 1) x L^{D) x L^{D) x 
L'^{D). We also need to estimate it using H2.11() - ()2.19() . Noticing that S* is independent 
of {9, q, T, S} and T* is only dependent on 9, we get 

S9t = 59yy -69- -f{y)[69 - 5T{y, 1)], < z < 1, (5.17) 

5qt = PrA5q 

- J((?\ V') + J{q^^p^)+PrRa{5Ty - 5Sy) + 59(1 - e^-')r],{z)), {y,z) € D, (5.18) 

5Tt = A6T-J{T^,iP^) + J{T^,^^) 

- J((l - e'-nrj,iz)9\^') + J((l - e'~nr^,{z)9\^P') - J((l - e'~^)7^,{z)So{y), 5iP) 

[((1 - e'-mi - 7(2/)) - e'-nVe{z) + (1 - e'^nVeiz) + 2e'-%{z) - e'-^v,]69 

- (1 - e'-')ve{zMy)6T{y, 1), (y, z) G D, (5.19) 
SSt = A6S - J{S\'^^) + J{S\ '^^) + J(S*, 5^), (y, z) G D. (5.20) 
The corresponding boundary conditions are: 

S9y{0,t) = 69yil,t) =0. (5.21) 

On the whole boundary: 

Sip = 0, 5Ai; = 6q = 0. (5.22) 
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At top z = 1: 

6T, + 6T\,=i = 0; 5S, = 0. (5.23) 

At bottom z = 0: 

6T^ = 5S, = 0. (5.24) 

At the lateral boundary y = 0,1: 

6Ty = 6Sy = 0. (5.25) 

The initial conditions are: 

600 = el - el 6qo = ql- ql 6To = - + (1 - e^-')ri,{z)60o, 5So = - Sl (5.26) 

where 56o = 56{y, 0), 6qo = 6q{y, z, 0),6Tq = 5T{y, z, 0) and 5Sq = 6S{y, z, 0). 
Multiplying H5.17() - (|5.2U|) by —5eyy, —A5q, —A6T and —A5S respectively, integrating over 
(0, 1) and D, we get 



1 d 
1 d 



■yi\^ = -\\S0yyf-\\69yf+ I \ {y)[6e - 6T {y , l)]5e yy dy , (5.27) 



2^I|V59|P = -Pr\\A6qf + J^iJ{q\^') - J{q^^'))A6q 

- PrRa [ {6Ty - 6Sy)A6q - / (1 - e^-')r],{z))6eA6q, (5.28) 
Jd Jd 

l±(\\V6Tf + \my,l)f) = -\\A5Tf+ [ {J{T\ij^) - J{TM^^))A5T 
2 dt Jd 

+ [ {J{{1 - e'-^)rj,{z)e\i;') - J{{1 - e'-^)ve{z)e\^'))A6T 



D 

[((1 - e'-')il - 7(y)) - e'-^Veiz) + (1 - e'-^rieiz) + 2e'~'r]'^iz) - e'-^r],]5eA6T 
+ / {l-e'-')r^,izhiy)6T{y,l)A6T 

JD 

+ f J((l - e^-')7],{z)So{y),64^)A6T, {y, z) G D, (5.29) 

JD 

^|||V5Sf = -||A<55f 
+ / {J{S\il:^)-J{S'^,il^^))A5S+ I J{S:,6iP)A5S, {y,z)eD. (5.30) 

JD JD 



Note that 
"1 







r{y)[5e - 6T{y, l)]Seyydy < -\\Seyyf + My)\\L^ fmf + \\dT{y, l)f ), 
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-PrRa / {6Ty - 6Sy)A6q < — \\Adqf + 2PrRa^{\\V6Tf + \\V6Sf ), 

/Pr 1 
^(1 - e'-^)rj,{z))6eA6q < —\\A6qf + —\m\ 

- I [((1 - - l{y)) - e'^nv.iz) + (1 - e'-^Vei^) + 2e'~^4i^) - e'-^ii,]5dA5T 

JD 

<\\\A5Tf + C\m\ 

/ (1 - e'-^)ii,{z)^{y)5T{y, l)A5T < hA5Tf + ^||5T(y, l)f, 
JD 6 2 

/ J((l - e'-^)v.iz)Soiy), 6^)A5T < hA6Tf + CXii\\Sof + \\S'£)\\qf, 
JD 6 

/ JiS:,6mSS < hA6Sf + CAidlFf + \\F'f)\\q\\\ 

JD 2 

Now by Lemma 13.11 we get 

/ iJiq\i;') - J{q\^P^))A5q < \\A6q\\\\Jiq\i;^) - Jiq^iJ^)\\ 
JD 

<^\\A6qf + Cil + X,n\75qf, 
[ {JiT\i;^) - J(r2, V'))A5T < \\A6T\\\\JiT\ij^) - J{T^ij')\\ 

JD 

<\\\A6Tf + Ci\\\75qf + \V6Tf), 
[ (J((l - e'-nv.{z)9\4'') - J((l - e'-^)rj,{z)9^i;^))A6T 

JD 

<\\\A6Tf + C{\\V6qf + \\V6ef), 
[ {JiS\^^) - J{S\^'^))A6S < l\\A5Sf + C{\\Vdqf + \V5Sf). 

JD I 

Putting the above estimtates together, we conclude 

< + ||5T(y, l)f + llV^gf + \\^bTf + (5.31) 

Taking N large enough, multiplying (|5.15p by and then adding to (|5.31() . we imply that 
there exists a positive constant a2 such that 

\j^mH? + 2NPrRa^\,{m? + \\bTf + \\bSf)) 

22 



+ llV^gf + \my, l)f + llV^Tf + \\V6S\\') 
< -a2{N\\6q\\^ + 2NPrRa^Xi{\\6e\\^ + \\6Tf + \\SSf)) 
- a2{\\Seyf + \\V6qf + \\V6Tf + ^\\V6Sf + \\5T{y, l)f). (5.32) 
By the Gronwall inequality, we have 

N\\6q\\^ + 2NPrRa^Xi{\\69\\^ + \\6Tf + \\SSf) 

+ \\S0yf + W^^Sgf + \\5T{y, l)f + \\V5Tf + \\S/5Sf 

< e-"2*(iV||(^gof + 2NPrRa^Xi{\\69o\\^ + ||<5Tof + ||(55of ) 

e-"'\\\Seoy\\^ + \\V6qof + \m{y, l)f + HV^Tof + ||V55of ). (5.33) 

This estimate tehs us that any two solution trajectories inside the absorbing set approach 
each other as time goes on. This is the so called strong contraction property. 

Remark 5.1 In fact, here we get the estimate is for |lV5T|p and ||V5S'|p = ||V5S'|p. But 
llVlrf = ||V((5r + (l-ei-^)??,(z)50)f , llV^f = llVfJSf and we have the estimate for 
\\60\\i. Hence we can obtain the similar estimate for the original 5T. 

Therefore we have the following theorem. 

Theorem 5.2 (^Strong contraction property^ Assume that the freshwater flux F{y) 
has zero mean as in \4-l^ ^''^d the ocean basin's latitudinal fraction function ^(y) is bounded 
as in \4-'^ - Let + H^alP + ||5'o||^2 + + ||-^||^2 be small enough. Then the coupled 
atmosphere- ocean system hl.l]) - ^T^\j has the strong contraction property. 

Now we come back to the issue of the time-almost periodic (in particular, time-periodic 
and time-quasi-periodic) motion in the coupled atmosphere-ocean system. First, we give 
the definitions about almost periodic function and pullback attractor. 
A function ip : M ^ X, where {X,dx) is a metric space, is called almost periodic ^ and 
if for every e > there exists a relatively dense subset of M such that 

dx {ip{t + T),^p{t)) < e 

for all t G iR and r € M^. A subset M Q M is called relatively dense in M if there exists 
a positive number / G JR such that for every a ^ ]R the interval \a,a + 1]{~\]R of length / 
contains an element of M, i.e. Mp|[a, a + /] 7^ for every a £ JR. 
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In order to study the temporally almost periodic solutions of H1.1I) ~ H1.9|) . we need some re- 
sults from the theory of nonautonomous dynamical systems. Consider first an autonomous 
dynamical system on a metric space P described by a group = {0t}teiR of mappings of 
P into itself. 

Let X be a complete metric space and consider a continuous mapping 

^: 1R+ X Px X ^ X 

satisfying the properties 

$(0,p, •) = idx, ^{T + t,p,x) = <!>{T,etp,<^{t,p,x)) 

for all t, T ^ JR^ , p & P and x € X. The mapping <I> is called a cocycle on X with respect 
to e on P. 

The appropriate concept of an attractor for a nonautonomous cocycle systems is the 
pullhack attractor. In contrast to autonomous attractors it consists of a family subsets of 
the original state space X that are indexed by the cocycle parameter set. 
A family A = {Ap}p^p of nonempty compact sets of X is called a pullback attractor of 
the cocycle $ on X with respect to 9t on P if it is (^-invariant, i.e. 

^{t,p,Ap) = Ae,p for all telR'^,peP, 

and pullback attracting, i.e. 

lim H*x e^tP, D),A„) = for all D G K{X), p £ P, 

where K{X) is the space of all nonempty compact subsets of the metric space {X,dx)- 

Here is the Hausdorff semi-metric between nonempty compact subsets of X, i.e. 

{{"^{A^B) := maXagA dist(a, -B) = maXa^A ^^'O-beB dx {A, b) for A, B £ K{X). 

The following theorem combines several known results. See Crauel and Flandoli [7j, Flan- 

doli and Schmalfufi ^1] , and Cheban |^ as well as [HI E] for this and various related 

proofs. 

Theorem 5.3 Let ^ he a continuous cocycle on a metric space X with respect to a group 
of continuous mappings on a metric space P. In addition, suppose that there is a 
nonempty compact subset B of X and that for every D G K{X) there exists a T[D) G 
, which is independent of p G P, such that 

<i>{t,p,D)cB for all t>T{D). (5.34) 

Then there exists a unique pullback attractor A = {Ap}p^p of the cocycle ^ on X , where 

Ap= [\ y ^{t,e_tp,B). (5.35) 
teiR+ 
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Moreover, if the cocycle $ is strongly contracting inside the absorbing set B. Then the 
pullback attractor consists of singleton valued sets, i.e. Ap = {a*{p)}, and the mapping p 
I— > a*{p) is continuous. 

The solution operators St.to for Hl.l() - ()1.9|1 form a cocycle mapping on X = H^{0,1) x 
Hq{D) X H^{D) X H^{D) with parameter set P = M, where p = to, the initial time, and 
dtto = to + t, the left shift by time t. However, the space P = M is not compact here. 
Though more complicated, it is more useful to consider P to be the closure of the subset 
{9tf, t € M}, i.e. the hull of /, in the metric space Ll^ {M, X)) of locally L2(iR)-functions 
f : ]R ^ X with the metric 



dp{f,g) ■■= f;2-^min|l,y^ 



N 

\\\fit)-g{t)Wdt 

N 



with 6t defined to be the left shift operator, i.e. 9tf{-) '■= /(• + t), where ||| • ||| denotes 
the norm in X. By a classical result ^ 02]) a function / in the above metric space is 
almost periodic if and only if the the hull of / is compact and minimal. Here minimal 
means nonempty, closed and invariant with respect to the autonomous dynamical system 
generated by the shift operators 9t such that with no proper subset has these properties. 
The cocycle mapping is defined to be the solution uj{t) = {9, q, T, S} of Hl.l|l - H1.9|) starting 
at loq = {6*0, qo,TQ, Sq} G X at time to = for a given forcing mapping f G P, i.e. 

^(i,/,wo) := loo, 

where we have included a superscript / on S to denote the dependence on the forcing 
term /. (This dependence is in fact continuous, see §3). The cocycle property here follows 
from the fact that S/i^wo = 'S'fl^p^o '^o for all t > to, to ^ M, uq £ X and f G P. 

Following Theorem 5.3 and the dissipativity and strong contraction properties shown in 
the last two sections, we conclude that the coupled atmosphere-ocean system Hl.l() - (|1.4|l 
has a unique pullback attractor, consists of the singleton valued component {a*{p)} G A 
and the mapping p a*{p) is continuous on P. As in Duan and Kloeden jJU] or Gao, Duan 
and Pu JZj, we now show that this singleton attractor a*{p) defines an almost periodic 
solution. 

In fact, the mapping p i-^- d*{p) is uniformly continuous on P because P is compact subset 
of Lf^^ {M, X)) due to the assumed almost periodicity. That is, for every e > there exists 
a 5(e) > such that ||a*(p) — a*(g)|| < e whenever dp{p,q) < 6. Now let the point p (= 
/, the given temporal forcing function) be almost periodic and for 6 = 6{e) > denote by 
Mg the relatively dense subset of M such that dp{9t+Tp, 9tp) < S for all r G Ms and t G 



25 



M. From this and the uniform continuity we have 



\\a*{et+rp)-a*i9tP)\\<e 

for all t e M and r G ^^(g). Hence t i-^ {6* , q* ,T* , S*}{t) := a*{6tp) is almost periodic, 
and it is a solution of the coupled atmosphere-ocean model. It is unique as the single- 
trajectory pullback attractor is the only trajectory that exists and is bounded for the 
entire time line. Therefore, we have the following result. 

Theorem 5.4 (T'eriodic, quasiperiodic and almost periodic motion^ Assume that 
the freshwater flux F{y) has zero mean as in i4-l\ ) ^-''^d the ocean basin's latitudinal frac- 
tion function ^{y) is hounded as in \4-^ - Let a} -|- HS'alP + ||<S'o|||^2 + + ||^|||^2 be 
small enough. Then the coupled atmosphere- ocean system il i. jjj - fTT^ has unique time- 
periodic, quasiperiodic and almost periodic motions, when the external fluctuation f in 
the atmospheric energy balance model is time-periodic, quasiperiodic and almost periodic, 
respectively. 

This result may be relevant to the El Nino-Southern Oscillation phenomenon. Ei Nino is 
a well-known climate phenomenon in the atmosphere-ocean system. Originally, it refers 
to a seasonal invasion, along the coast of Peru around Christmas, of a warm southward 
ocean current that displaced the north-flowing cold current. Today, it is regarded as a part 
of a phenomenon called El Nino-Southern Oscillation (ENSO), a continual but "quasi" - 
periodic, perhaps irregular, cycle of shifts in ocean and atmosphere condition that affect 
the globe climate PHIHIOH]. 

6 Summary 

We have investigated the dynamical behavior of a coupled atmosphere-ocean system. 
First, we show that the coupled atmosphere-ocean system is stable under the external fluc- 
tuation in the atmospheric energy balance relation (Theorem 13. 2() . Then, we estimate the 
atmospheric temperature feedback in terms of the freshwater flux, heat flux and the exter- 
nal fluctuation at the air-sea interface, as well as the earth's longwave radiation and the 
shortwave solar radiation (Theorem 14.1)1 . Finally, we prove that the coupled atmosphere- 
ocean system has time-periodic, quasiperiodic and almost periodic motions (under suitable 
conditions on the physical quantities such as freshwater flux, the earth's longwave radiative 
cooling coefficient and the shortwave solar radiation profile), when the external fluctua- 
tion in the atmospheric energy balance relation is time-periodic, quasiperiodic and almost 
periodic, respectively (Theorem 15. 4j) . 
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